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ABSTRACT 

In order to assess the contribution of Lyman break galaxies (LBGs) and Ly- 
man alpha emitters (LAEs) at rcdshifts 3 < z < 7 to the ionization of intergalactic 
medium (IGM), we investigate the escape fractions of ionizing photons from supernova- 
dominated primordial galaxies by solving the three-dimensional radiative transfer. The 
model galaxy is employed from an ultra-high-resolution chemodynamic simulation of 
a primordial galaxy by Mori & Umemura (2006), which well reproduces the observed 
properties of LAEs and LBGs. The total mass of model galaxy is 10 11 M Q . We solve 
not only photo-ionization but also collisional ionization by shocks. In addition, accord- 
ing to the chemical enrichment, we incorporate the effect of dust extinction, taking 
the size distributions of dust into account. As a result, we find that dust extinction 
reduces the escape fractions by a factor 1.5 — 8.5 in the LAE phase and by a factor 
2.5 — 11 in the LBG phase, while the collisional ionization by shocks increases the 
escape fractions by a factor w 2. The resultant escape fractions are 0.07 — 0.47 in 
the LAE phase and 0.06 — 0.17 in the LBG phase. These results are well concordant 
with the recent estimations derived from the flux ratio at 1500 A to 900 A of LAEs 
and LBGs. Combining the resultant escape fractions with the luminosity functions of 
LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z = 3 — 5. 
However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM 
at z > 6. That implies that additional ionization sources may required at z > 6. 

Key words: radiative transfer - ISM: dust, extinction - galaxies: evolution - galaxies: 
formation - galaxies: high-rcdshift 



1 INTRODUCTION 

One of momentous issues regarding the evolution of inter- 
galactic medium (IGM) is the ionization history of the uni- 
verse , which significantly influences the galaxy formation 
(e.g.,ISusa fc Umemural200d;IUmemura. Nakamoto. fc Susal 
l200ll : ISusa fc Umemurall2004l ). The Wilkinson Microwave 
Anisotropy Probe (WMAP) results provi de a wealth of in- 
formation about the cosmic reionization (|Page et alj|2007t 
iDunklev et alj|2009l ). However, the detailed history of reion- 
ization and the nature of ioniz i ng so urces are not yet fully 
understood. lHaardt fc Madaul (| 19961 ) pointed out that the 
UV backg r ound radiation is dominated by quasars at z < 4. 
iFan et alj (|200ll ) showed, using the SDSS sample, that the 
bright-end slope of the quasar luminosity function at z >, 4 
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are considerably shallower than that at low-redshifts, and 
they concluded that quasars cannot maintain the ioniza- 
tion of IGM at z >j 4. Subsequently, a lot of arguments 
have been concentrated on the possibility that the IGM is 
ionized mainly by UV radiation from high-z star-forming 
galaxies like Lyman break galaxies (LBGs) and Lyman a 
emitters (LAE s) (e.g.. lFan et alj|200d : feouwens et al.ll2007l ; 
lGnedidl200sl ). However, the estimate of the contribution 
of LBGs or LAEs suffers significantly from the ambigu- 
ity regarding the escap e fractions of ionizing photon s from 
star-forming galaxies ([Razoumov fc Sommer-Larsenl 120061 ; 
iGnedin, Kravtsov fc Chenll2008l ). 

Observationally, the escape fractions of ionizing photons 
are assessed by the relative escape fractions /esc,rei> using the 
flux ratio at 1500 A to 900 A, (F1500/-F900) obs , which are 
defined by 



y*esc,rel 



(L1500/L900)mt 
(F1500/F900)obs 



exp(r 900 ) 



(1) 
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where TgJg> M represents the line-of-sight opacity of the IGM 
for 900 A photons. Normally, the intrinsic luminosity ratio 
at 1500 A to 900 A, (L1500/L900) in t, is assumed to be 3 as 
a fiducial value. ISteidel. Pettini fc Adelbergerl (|200ll ) found 
/csc.roi Si 0. 5 from the composite s pectr um of 29 LBGs at 
z ~ 3, and iGiallongo et all (|2002T ) and llnoue et all (|2005h 
estimated the upper limit of / CS c,rci 0.1 — 0.4 for some 
LBGs at z ~ 3. The direct detection of ionizing photons 
from high-z star-forming galaxies has been accomplished re- 
cently as a consequen ce of intensive and continuous exertion. 
IShaplev et al.l (|2006l ) detected the escaping ionizing photons 
from 2 LBGs in the SSA22 field at z = 3.1 and they esti- 
mated the average rel a tive e scape fraction / CS c,roi = 0.14. 
Moreover, llwata et al.l (|2009l ) successfully detected the Ly- 
man continuum emission from 10 LAEs and 7 LBGs within 
197 samples of LAEs and LBGs in the SSA22 field. They 
have shown that the mean value of relative escape fractions 
for 7 LBGs is 0.11 after a correction for dust extinction, and 
can be 0.20 if IGM extinction is taken into account. 

Theoretically, the accurate estimation of escape frac- 
tions requires the three-dimensional (3D) radiative transfer 
calculations of ionizing photons traveling in inhomogeneous 
interstellar medium. Such 3D radiative transfer calculations 
are fairly recent attempts, beca use a great deal of com- 
putat ions are demanded. First, ICiardi. Bianchi fc Ferraral 
|2002r ) calculated the escape fractions in simple model 
clouds, where the density distributions are smoothed Gaus- 
sian or fractally inhomogeneous. More recently, using gas- 
dynamical simulations of galaxy formation, the escape 
fractions are estimated by s o lving 3D r adiative transfer. 
iRazoumov fc Sommer- Larsenl l|2006l . l2007h have found that 
the escape fractions decline from several per cent at z = 3.6 
to 0.01-0.02 at z = 2.39, due to higher gas clumping at 
lower redshifts. But, in these simulatio ns, the effects o f dust 
extinction are not taken into account. iGnedin et al.l (|2008l ) 
performed a cosmological simulation on the formation of a 
disk-like galaxy that provides the three-dimensional distri- 
butions of absorbing gas in and around the galaxy. Then, 
3D radiative transfer was solved including the dust extinc- 
tion. As a result, the escape fractions turned out to be as 
low as a few per cent. In this model, the bulk of stars is 
embedded deep inside the optically-thick HI disk and there- 
fore most of ionizing photons emitted from hot stars are 
absorbed by neutral hydrogen in the HI disk. Only a small 
fraction of ionizing photons from the stars that are located 
near the edge of disk can escape from the galaxy. Hence, 
almost regardless of the effects of dust extinction, the re- 
sultant escape fractions become quite small. It implies that 
star-forming galaxies cannot give a significant contribution 
to the IGM ionization at z > 3. For higher redshift z > 8, 
on the other hand, some numerical simulations have already 
shown that large escape fractions of some tens per cent 
are p ossible for low-mass galaxies, which include Pop HI 
stars dAlvarez. Bromm fc Shapirdl200d;lwhalen et al.ll2004l : 
iKitavama et al.ll2004l : IWise fc Cenll2009l 'l . However, in order 
to assess the contribution of LAEs and LBGs observed at 
3 <, z < 7 to IGM reionization, we should evaluate the UV 
escape fractions from more massive galaxies. The escape of 
UV photons can sensitively depend on the gravitational po- 
tential (Whalen et al. 2004, Kitayama et al. 2004). Also, dust 
extinction should be taken into account in such primordial 
galaxies. 



In this paper, we reconsider the escape fractions of ion- 
izing photons from high-z primordial galaxies by employ- 
ing a su pernova-dominated primo rdial galaxy model pro- 
posed by iMori fc Umemural (|2006l ) , which well reproduces 
the observed properties of LAEs and LBGs. The total mass 
is lO n M0, and the star formation with the Salpeter IMF is 
included. In this model, stars are distributed more extend- 
edly in the galaxy and the bulk of interstellar gas is collision- 
ally ionized by supernova-driven shock. The spread of stellar 
distributions and the shock heating can lead to diminishing 
the absorption of ionizing photons and may enhance escape 
fractions. Here, we perform the 3D radiative transfer cal- 
culations including not only the photo-heating but also the 
shock heating. We incorporate the dust extinction that is 
consistent with the chemical evolution of galaxy. Then, us- 
ing the resultant escape fractions, we explore whether LAEs 
and LBGs can contribute to the ionization of IGM in a high- 
z universe. In the present analyses, the cosmological param- 
eters are assumed to be Ho = 70 km s _1 Mpc -1 , Qm = 0.3 
and S1a = 0.7. In §2, the model and numerical method are 
described. In §3, the numerical results on escape fractions 
are provided. In §4, we argue the contributions of LAEs and 
LBGs to the IGM ionization. In §5 is devoted to the sum- 
mary. 



2 MODEL AND METHOD 

As a model galaxy, we adopt the high-resolution hy- 
drodynamic simulations (1024 3 fixed Cartesian grids) by 
IMori fc Umemural (|2006l 'l. which are coupled with the col- 
lisionless dynamics for dark matter particles as well as star 
particles and also the chemical enrichme nt. The simulation 
box s ize is 134 kpc in physical scales. I Mori fc Umemural 
l|200q ) demonstrated that an early proto-galactic evolution 
with multitudinous type II SNe (SNell) exhibits intense Ly- 
man q emission, well resembling LAEs. Subsequently, the 
galaxy shifts to a stellar continuum radiation-dominated 
phase, which appears like LBGs. In the present analysis, 
the stages from t agc = O.lGyr to f ago = 0.3 Gyr are defined 
to be the LAE phase, and the stages from t agc = 0.5 Gyr to 
£a g c = 1.0 Gyr are the LBG phase. 

As a post process, we calculate the escape fractions of 
ionizing photons with an accurate radiation transfer scheme. 
The data of the hydrodynamic simulations are coarse- 
grained into 128 3 Cartesian grids to solve radiation trans- 
fer. We use t he Authentic Radiation Transfer meth od (ART) 
developed bv lNakamoto. Umemura fc Susa (|200ll ). The per- 
formance of this scheme has already report ed as a part of 
the comparison study bv llliev et al. (|2006al ). The radiation 
transfer equation is solved along 128 2 rays with uniform an- 
gular resolution from each source. At a point of optical depth 
r, the specific intensity is given by I v (t) = I V {Q) exp (— r), 
where 7„(0) is the intrinsic intensity of ionizing radiation 
and r is the optical depth of neutral hydrogen and dust. As 
for scatt ering photons, we employ the on-the-spot approxi- 
mation (|Osterbrocklll989T ). in which scattering photons are 
assumed to be absorbed immediately on the spot. We ob- 
tain the ionization structure assuming the ionization equi- 
librium, r 7 7iHi + r C 7iHi«c = ceBn p n c , where r 7 ,r c and 
«b are the photo-ionization rate, the collisional ionization 
rate and the recombination rate to all excited states, re- 
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spectively. We continue the radiative transfer calculation 
recurrently until the ionization structure converges. In a 
supernova-dominated model galaxy, the collisional ioniza- 
tion occurs mostly in low density, high temperature regions 
with n < 10 -3 cm -3 and T co ii > 10 4 K. In such regions, stel- 
lar UV radiation contributes to the increase of ionization 
degree, but does not much to the increase of the temper- 
ature. On the other hand, neu tral regions are photoion- 
ized and heated up to ~ 10 4 K (|Umemura fc Ikeuchi|[l98i ; 



iThoul fc Weinberg! Il996h . when ionizing radiation is irra- 
diated. Hence, in the present analysis, we assume the gas 
temperature to be T = max{T co n, 10 4 K} in ionized regions, 
and we do not update the temperature in iteration of ra- 
diation transfer calculation. To evaluate the escape fraction 
of ionizing photons, we count all photons above the Lyman 
limit that escape from the calculation box. 

The number of ionizing photons emitted from source 
stars is computed based on the theoretical spectral en- 
ergy distribution (SED) give n by a population syn- 
thesi s scheme, PEGASE v2.0 H ioc fc Rocca-Volmerangg 
1997). We assume the Salpeter (1955) initial mass func- 
tion in the mass range of Q.IMq — 50Mq. As for dust 
grains, we adopt the empirical size dist ribution nd(ad) oc 
a~ 3 - 5 (|Mathis. Rumpl fc Nordsieck|[i977h in the range from 
0.1/im-1.0^m , where ad is the radius of a dust grain. We 
assume refractory grains like silicates, for which the pho- 
todestruction of dust by UV radiation is neg ligible over 
the Hubble time-scale (jDraine fc Salpeterlll979h . The dust 
grains are distributed proportionally to the metallicity cal- 
culated in the hydrodynamic simulations with the relation 
of md = 0.01m 9 (Z/Zq), where mj, m g , and Z are the dust 
mass, gas mass, and metallicity in a grid. The density in a 
dust grain is assumed to be 3 g cm -3 like silicates. The dust 
opacity is given by drd us t = Q(v)na\nads, where Q(v),a,d 
and rid are the absorption Q-value, dust size and number 
density of dust grains, respectively. Since the assumed range 
of dust size is larger than the wavelengt h of Lyman limit, w e 
assume Q(v) = 1 for ionizing photons (|Draine fc Lee|[T984h . 



3 RESULTS 

3.1 Galactic Evolution 

Figure 1 shows the snapshots for the evolution of model 
galaxy as a function of redshift from t agc = 0.1 to 1.0 Gyr. 
The density of gaseous a nd st ellar components are adopted 
from iMori fc Umemural (|200rJ l . The dust density is evalu- 
ated by the metallicity. As a result of multiple supernova 
explosions, dust is distributed more extendedly than stars, 
which is significantly relevant to the absorption of ioniz- 
ing photons. Bottom panels show the calculated ionization 
structure in terms of neutral hydrogen fractions xhi- 

First, we see the evolution of the model galaxy. At 
tagc =0.1 Gyr, stars form in high-density peaks in sub- 
galactic condensations and the burst of star formation starts. 
Then, massive stars in the star-forming regions explode as 
SNell one after another. The gas in the vicinity of SNell 
is quickly enriched with the ejected heavy elements and 
then interstellar dust is locally procreated. However, a large 
amount of gas still retains a metal-free and dust-free state. 
The spatial distributions of heavy elements are highly inho- 
mogeneous, where gas enriched with —5.1 <, [Si/H] <, — 1.1 



and —5 <, [O/H] <, — 1.0 coexist with virtually metal-free 
gas. Supernova-driven shocks collide with each other to gen- 
erate large-scale hot 10 6 K) bubbles reaching a higher 
ionization degree. At t ag e = 0.3 Gyr, roughly 50 per cent of 



the total volume is highly ionized with xhi = 10 



10" 



where the ionization degree is controlled by the collisional 
ionization by shock heating and the photo-ionization by UV 
radiation from hot young stars. The spatial distributions of 
the heavy elements are still highly inhomogeneous in the 
range of -2.4 < [Si/H] < - 0.4 and -2.5 < [O/H] < - 0.5. 

After t ag c = 0.5 Gyr, the hot bubbles expand and 
blow out into the intergalactic space. As a result, more 
than 80 per cent of volume is occupied with highly ion- 
ized gas with xhi = 10~ 6 ~ 10~ 7 . Finally, the merger of 
sub-galactic condensations promotes the mixing of heavy- 
elements and weakens the spatial inhomogeneities of heavy- 
element abundance and ionization degree. As a result, the 
heavy-element abundance of interstellar medium (ISM) con- 
verges to -0.4 < [Si/H] < 0.1 and -0.3 < [O/H] < 0.2 with 
small dispersion, and eventually 95 per cent of volume is 
filled with ionized gas. This means that the metallicity 
reaches around the solar abundance in 10 9 yr. It is consistent 
with the previous works on the ellipt i cal galaxy formation 
llArimoto fc Yoshiilll987|; iGibsonll 19971 ; iKodama fc Arimotol 
ll997tlKawata fc Gibson! l2003al [bj). The mean value of mass 
weighted heavy-element abundance is [Si/H] ~ —1.0 and 
[O/H] ~ -0.9 at the LAE phase, while [Si/H] ~ -0.4 and 
[O/H] ~ —0.3 at the LBG phase. These are translated into 
the mass weighted metallicity as 0.14 Zq at LAE phase and 
0.52 Zq at LBG phase, which a re co ncordant with th e obser - 
vations bv lPettini et al.l (120011 ) and iMannucci et al.l (|2009h . 

The luminosity of Lyman a emission, which is the cool- 
ing radiation by interstellar gas, reaches 2.0 x 10 43 erg s _1 
at t ag c = 0.1 Gyr and 1.6 x 10 43 erg s _1 at t agc = 0.3 Gyr, 
respectively. They nicely match the observed luminosity of 
LAEs and also well resemble LAEs with respect to other 
properties. After t ag e ^ 0.5 Gyr, the Lyman a luminosity 
quickly declines to several 10 41 erg s -1 that is lower than the 
detection limit. Then, the galaxy shifts to a stellar contin- 
uum ^adjatioir^d^miMtecFphase, which appears like LBGs 
(see IMori fc Umemurall2006l ) . 



3.2 Escape Fractions 

Figure 2 shows the time evolution of escape fractions of ion- 
izing photons derived by the full radiation transfer calcula- 
tions. The upper panel represents the absolute escape frac- 
tions, which are defined by 



/- 



(2) 



total 

7 



where N^ otal is the total number of pristine photons radi- 
ated from stars and N2, c is the number of photons escaped 
from the simulation box. Filled circles in the upper panel 
of Fig. 2 represent the resultant escape fractions / eS c with 
dust extinction. In the LAE phase, / csc = 0.47 at t agc = 0.1 



Gyr and f C! 



0.23 at t a 



0.3 Gyr, while in the LBG 



phase, / osc = 0.19 at t agc = 0.5 Gyr and / OS c = 0.17 at 
tagc = 1.0 Gyr. The escape fractions have the dependence 
on the size distributions of dust, since smaller grains result 
in lager extinction for a given dust-to-gas ratio. If we ex- 
tend the range to smaller grains as 0.02 — l.O^im, the escape 
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Figure 1. Snapshots of the evolution of the model galaxy with total mass of 1O 1:L M0, at t age = 0.1,0.3,0.5 and 1.0 Gyr. Each panel in 
row corresponds the spatial distributions of gas density p s (g cm -3 ), dust density (g cm -3 ), stellar density p s (g cm -3 ), and fractions 
of neutral hydrogen in logarithmic scale Xhi> respectively. The simulation box is 134 kpc in physical scale. 



fractions are reduced to / csc = 0.1 — 0.3 in the LAE phase 
and / C sc = 0.09 — 0.1 in the LBG phase. If the grains are 
even smaller as 0.03 — 0.3/im, the escape fractions become 
/ CS c = 0.07-0.22 in the LAE phase and / e8C = 0.06-0.07 in 
the LBG phase. Thes e values are concord ant w ith the esti- 
matio ns for LBGs bv lShaplev et all (|2006t ) and llwata et alj 
l|2009D . 

The escape fractions are significantly regulated by in- 
terstellar dust. In Fig. 2, we also show the dust-free model 
(open circles), where the absorption by interstellar dust is 
artificially switched off. We find that the dust extinction re- 
duces /esc by a factor of 1.5-8.5 in the LAE phase. On the 
other hand, The reduction by dust extinction is a factor of 
2.5 — 11 in the LBG phase. The change of the reduction 
factor is linked to the enrichment of heavy elements in the 
galaxy. 



iGnedin et ah | (|2008l) have shown that the escape frac- 
tions are as small as a few per cent, which is much smaller 
than t he present results. In the calculations bv lGnedin et al.l 
(|2008l ), most of young stars are distributed in the dense re- 
gion of galactic disk. In our model, stars are extendedly dis- 
tributed, and also the bulk of interstellar gas is collision- 
ally ionized by the SN shock heating. Since the interstellar 
medium is moderately optically-thin for ionizing photons, 
ionizing photons can escape through collisionally-ionized re- 
gions. To assess the effects of collisional ionization, we have 
tentatively calculated a case assuming T = 10 4 K and no 
collisional ionization. As a result, we have found that the 
escape fractions are reduced to be / CS c = 0.21 at t aga = 0.1 
Gyr and / CS c = 0.13 at t agc = 1.0 Gyr. This implies that the 
collisional ionization by shocks can contribute to enhance 
the escape fractions by a factor of two. 



Escape of ionizing photons 5 



0.0 
1.0 





Gyr 

Figure 2. Upper panel : Absolute escape fractions / eS c for the 
simulated galaxy as a function of evolutionary time. Blue symbols 
represent the escape fractions / csc with dust extinction, where 
blue filled circles, open diamonds, and open pentagons show the 
escape fraction for dust size of 0.1 — 1.0/xm, 0.02 — l.O^tm, and 
0.03 — 0.3/jm, respectively. Black open circles denote / csc with- 
out dust extinction. The black open triangles show the dust sput- 
tering model. Lower panel : Absolute escape fractions / esc for 
dust size of 0.1 — 1.0pm are compared to relative escape fractions 
/csc rel that are derived from the flux ratio at 1500 A to 900 A 
[see equation JTJ]. The filled circles represent the absolute escape 
fractions / eS c and the open squares do the relative escape frac- 
tions /csc.rcl- The variations by viewing angles are represented by 
a pink belt with the standard deviation shown by a narrower red 
belt. 



If dust grains are in high-temperature regions with 
T > 10 6 K, they may be destructed by sputtering pro- 
cess. The grain radius ad decre ases at a rate dad/ dt 

~3Q 



\2.5l 



dDraine fc Salpeterl fl979l : 



iTsai fc Mathews! Il995l ; iMathews fc Brighentill2003ri . where 
/i w = 3.2 x 10" 18 cm 4 s" 1 and T s = 2 x TO 6 K are 
fitting parameters, and n p is the proton number density. 
Therefore, the destruction time scale of ~ 0.1 /j, m dust is 
shorter than the time scale of galaxy evolution ( < lGyr), 
if n p >, 10~ 4 cm~ 3 and T > 10 6 K. Since the present model 
galaxy has n p > 10~ 4 on average in high-temperature re- 
gions, we roughly suppose that all dust in high-temperature 
regions (T > lO 6 ^) is evaporated by sputtering. Open tri- 
angles in the upper panel of Fig. 2 show the escape frac- 
tions when taking the sputtering into account. After includ- 
ing dust sputtering, / csc = 0.58 at t agc = 0.1 Gyr, and 
/csc = 0.17 at t ag c = 1-0 Gyr. At higher redshifts, the dust 
sputtering somewhat enlarges the escape fractions, because 
a part of interstellar dust is distributed in high-temperature 
regions. On the other hand, at lower redshifts, the results 



are basically the same as no-sputtering cases, since ionizing 
photons are mainly absorbed by dust in low-temperature 
star-forming regions. Although the sputtering model here 
may overestimate the destruction of dust, the escape frac- 
tions do not change very much by including the dust sput- 
tering. Also, in our simulations, the dust temperature does 
not rise over 100 K by photo-heating. Since the evaporation 
temperature of silicate dust is ~ 10 3 K, the photodestruction 
of dust is unimportant. 

In order to evaluate the relative escape fractions de- 
fined by (UJ), we make a mock observation of the simulated 
galaxy using the flux ratio at 1500 A to 900 A. Assuming 
(L1500/L900) int = 3 (|Steidel et al.ll200lh and Q(1500 A) = 
1, and also that the radiation flux at 1500 A is absorbed 
only by interstellar dust, we compute the transport of the 
radiation fluxes at 900 A and 1500 A. The resultant relative 
escape fractions /csc.rci are shown by open squares in the 
lower panel of Fig. 2. It turns out that the predicted rela- 
tive escape fractions (/ eS c,rci ~ 0.1 — 0.3) mat ch the average 
of ob served relative escape fractions as 0.14 (|Shaplev et al.l 
l2006h . Al so, it is consis t ent w it h the observational esti mates 
given bvlSteidel et all ri200ll ), iGiallongo et~aLl (|2002T l. and 
llnoue et all (|2005l . l200fj ). 



4 DISCUSSION 



IShaplev et alj (|2006l 'l estimated /csc.rci for 14 LBGs in the 
SSA22a field, and they detected ionizing photons from only 
2 LBGs, that is, C49 and D3. The reported relative escape 
fractions of 2 LBGs are extremely high as /esc,rei(C49) = 
0.65, and / C sc,rci(D3) ^ 1. 0. For the o t her 1 2 LBGs, only 
upper limits are suggested. Ilwata et al. | l|2009h reported the 
results of Subaru/Suprime-Cam deep imaging observations 
of the same field. They detected ionizing radiation from 7 
LBGs as well as from 10 LAE candidates. They also showed 
the large scatter of the observed UV to Lyman continuum 
flux density ratios. For the seven detected LBGs, the ratio 
ranges from 2.4 to 23.8 with a median value of 6.6. Then, the 
relative escape fractions for 7 LBGs ranges from 0.03 to 0.30 
after a correction for dust extinction, and can be 0.05 — 0.55 
if IGM extinction is taken into account. In addition, some 
of the detected galaxies show significant spatial offsets of 
ionizing radiation from non-ionizing UV emission. 

The spatial distributions of interstellar dust are 
highly inhomogeneous in LAEs, because it should follow 
the inhomogeneous heavy-element distributions (see also 



Mori. Ferrara fc Madaul |2002| ; |M ori. Umemura fc Ferraral 



2004h . Therefore, the observed escape fractions strongly de- 



pend on viewing angles. We compute the probability dis- 
tribution function, dp(f)/df = A r (/)/7V tota i, as a func- 
tion of relative escape fraction along a line of sight, where 
A^totai = 128 2 is the total number of bins in viewing angles 
and N(f) is the number of angular bins with a given escape 
fraction /. Figure 3 shows the probability distribution func- 
tion from redshift t ago = 0.1 to 1.0 Gyr. This figure clearly 
shows that the relative escape fractions vary by viewing an- 
gles significantly. For example, depending on viewing angles, 
the relative escape fractions / eS c,rei vary from 0.18 to 0.43 
at iage = 0.1 Gyr and from 0.12 to 0.47 at t agc = 1.0 Gyr. 
These probability distributions are indicated by a pink belt 
in the lower panel of Fig. 2. The standard deviation at each 
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Figure 3. Probability distribution of relative escape fractions. 
Each line corresponds to the probability distribution at t ago = 
O.lGyr (green), 0.3 Gyr (blue), 0.3Gyr (magenta) and 1.0 Gyr 
(red), respectively. 



redshift is also shown by a red belt. This effect allows some 
of observed L AEs and LBGs to h ave as high escape fraction 
as reported bv llwata et al.l l|2009l ). 

The present analyses show the absolute escape frac- 
tions /esc for LAEs and LBGs can be as large as / osc >, 0.17. 
Hence, LAEs and LBGs are potential sources for the IGM 
ionization at z >, 4. Here, we quantify the contributions of 
LAEs and LBGs to the IGM ionization. Specifically, we as- 
sess the emission rate of ionizing photons from LAEs and 
LBGs per unit comoving volume. The emission rate is eval- 
uated from the star formation rate based on the observed 
luminosity functions, with coupling the escape fractions for 
LAE and LBG phases obtained in the present analysis. We 
use the average escape fraction < / csc >= 0.35 for the LAE 
phase and < / csc >= 0.18 for the LBG phase. 

In Figure 4, the emission rate of ionizing photons 
per comoving Mpc 3 is shown as a function of redshift. 
Open symbols depict the emission rate for the samples 
of LAEs that listed in the caption. Filled symbols show 
the emission rate for LBGs, which is estimated by ex- 
trapolating the lum inosity function to L = 0.1L*_ 3 from 
ISteidel et all (Il999h with dus t extinction of E(B-V)=0.13 
(see also lYoshida et al.1 12006T). As a theoretical criterion, 
we adopt that by Madau. Haadt fc Reesl f|l999l ') . where the 
emission rate of ionizing photons required to balance the re- 



combination is given by N ion = 10 47 ' 4 C(1 + z) 3 s _1 Mpc" 
C is the clumping factor which param eterizes the inhomo - 
geneity of ionized hydrogen in the IGM. iMadau et all (jl999T ) 
adopted C = 30, based on the value co mputed by a cosmo- 
logica l radiative transfer simulation of iGnedin fc Ostrikerl 
(|l997l ). l0uchi et all (|2004l ) also adopted C = 30 when com- 
puting the number of ionizing photons needed to keep the 
IGM highly ionized at z = 5, as do a number of other au- 
thors. Although C may have some uncertainty, we assume 
C = 30 here as a fiducial value. 

Figure 4 clearly illustrates that observed LBGs can pro- 
vide the majority of ionizing photons at z = 3 — 5, and 
play an important role to keep the universe ionized. How- 
ever, LAEs are not capable of ionizing large volumes at the 
redshift z — 3 — 5. On the other hand, ionizing photons 
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Figure 4. The evolution of emission rate of ionizing photons 
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0.35 which is mean escape fraction at LAE phase. The blue 
filled symbol s repr esent the Nj on of LBGs derive d from lb) 
ISteidel et"al] dl999l). 2b) lYoshida et al.1 ll2006h. 3b) l lwata et all 



j2003l). 4b) iBouwens et al.N2006l) . 5 b) lOuchi et al.1 |2004l), 6b) 
ISawicki fc Thompson! l|2006l ). and 7b) iGabasch et al.l \20o4 ) with 
< /csc >= 0.18 which is mean escape fraction at LBG phase. 
The horizontal and vertical error-bars are arisen from the uncer- 
tainty of observations and the variation of escape fractions (LAE : 
/ osc = 0.22-0.47, LBG : / csc = 0.17-0.19), respectively. A solid 
line and a dotted line indicate the emission rate re quired to ion- 
ize th e IGM with C = 30 and C = 1, respectively ijMadau et al.l 
|l999f ). A dashed line represents the emission rate evaluated by 
the QSO luminosity function. 



from not only LBGs but also LAEs are not enough to ion- 
ize the IGM at z >, 6. Most of photons that ionize the uni- 
verse may come from u ndetected faint LAEs and LB Gs or 
other sources. Recently, IChoudhurv fc Ferraral (|2007l ) stud- 
ied the cosmic reionization history to account for a number 
of observational data, and pointed out that low-mass galax- 
ies hosting Pop III stars can be predominant ionizing sources 
of the IGM at high-z. Our results advocate their model. In 
the present analysis, all Lyman a photons are assumed to 
escape. Hence, the contribution of LAEs may be underesti- 
mated. In the future work, we intend to include Lyman a 
line transfer to compare the numerical results more precisely 
with the observations. 



5 SUMMARY 

We have performed three-dimensional radiation transfer cal- 
culations, based on a high-resolution hydrodynamic simula- 
tion of a supernova-dominated primordial galaxy, to obtain 
the ionization structure and explore the escape fractions of 
ionizing photons from LAEs and LBGs at high redshifts. 
The effect of dust extinction is incorporated according to 
the chemical enrichment, taking the size distributions of 
dust into account. As a result, we find that dust extinc- 
tion reduces the escape fractions by a factor of 1.5 — 8.5 
in the LAE phase and by a factor of 2.5 — 11 in the LBG 
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phase. The resultant escape fractions are 0.07 — 0.47 in the 
LAE phase and 0.06 — 0.17 in the LBG phase. These results 
are well concordant with recent observations. We have found 
that the combination of diffuse distributions of stars and su- 
pernova shock-heating is important for UV escape fractions. 
In the present galaxy model, young stars are extendedly dis- 
tributed, and also the bulk of interstellar gas is collisionally 
ionized by supernova shock-heating. Since the interstellar 
medium is moderately optically-thin and quite bubbly, ioniz- 
ing photons can escape through collisionally-ionized regions. 
The collisional ionization by shocks contributes by a factor 
of ~ 2 to the increase of the escape fractions. 

The relative escape fractions derived by mock observa- 
tions of the simulated galaxy match quite well the estimates 
by recent observations for LAEs and LBGs. To assess the 
contribution of LAEs and LBGs to the IGM ionization, the 
resultant escape fractions have been combined with the lu- 
minosity functions of LAEs and LBGs. As a result, we find 
that high-z LAEs and LBGs can ionize the IGM at z = 3 — 5. 
However, ionizing radiation from LAEs as well as LBGs is 
not enough to ionize the IGM at z > 6. That implies that 
undetected faint LAEs and LBGs or additional ionization 
sources may deter mine the IGM ioniz ation at z > 6. 

Very recently, IWise fc Cenl \200$ ) performed 3D radia- 
tion hydrodynamic simulations to assess the contribution of 
dwarf galaxies to cosmic reionization at redshift z = 8. They 
studied the UV escape fractions for low-mass galaxies in the 
mass range of Aftotai = 3 x 10 6 - 3 x 10 9 M G . As a result, they 
have shown that the UV escape fractions can reach up to ~ 
0.8 without dust extinction in halos with > lO 8 A/0 for a 
top-heavy initial mass function (IMF). However, in order to 
assess the contribution of LAEs and LBGs to IGM reioniza- 
tion, we should evaluate the UV escape fractions from more 
massive galaxies with dust extinction. Here, we have shown 
that a high-mass, metal-enriched galaxy at low redshifts can 
allow escape fractions as large as some tens per cent . We 
have found that an order-of-magnitude increase in the es- 
cape fractions can be attributed to the diffuse distribution 
of stars, while the collisional ionization further raises them 
by some factor. 
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